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Abstract 


The first study of W boson production in pPb collisions is presented, for bosons de¬ 
caying to a muon or electron, and a neutrino. The measurements are based on a data 
sample corresponding to an integrated luminosity of 34.6 nb^^ at a nucleon-nucleon 
centre-of-mass energy of = 5.02 TeV, collected by the CMS experiment. The 

W boson differential cross sections, lepton charge asymmetry, and forward-backward 
asymmetries are measured for leptons of transverse momentum exceeding 25 GeV/c, 
and as a function of the lepton pseudorapidity in the |?/iab| < 2.4 range. Deviations 
from the expectations based on currently available parton distribution functions are 
observed, showing the need for including W boson data in nuclear parton distribu¬ 
tion global fits. 
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1 Introduction 

Electroweak boson production in proton-nucleus and nucleus-nucleus collisions at the CERN 
LHC offers a unique opportunity to probe nuclear parton distribution functions (nPDEs) 
Leptonic decays of electroweak bosons are of particular interest since leptons do not inter¬ 
act strongly with the medium produced in these collisions JSlIS]. As compared to those in a 
proton, the nPDEs are expected to be depleted (shadowing) for partons carrying small momen¬ 
tum fractions x < 10^^, and enhanced (anti-shadowing) in the 5 x 10^^ ^ ^ ^ 10^^ range 17]. 
However, because of the lack of available data, parton densities are less precisely known for 
nuclei than for nucleons. As a consequence, precise calculations describing hard processes in 
high-energy heavy ion collisions are limited by uncertainties in the nPDEs. Eor W boson pro¬ 
duction, the dominant processes at LHC energies are ud —> W+ and du —> W^, principally 
reflecting interactions that take place between valence quarks and sea antiquarks. According 
to Ref. |4|, PDE nuclear modifications could affect the yield of W bosons in pPb collisions at 
the LHC by as much as 15% in certain kinematic regions. Therefore, precise measurements 
of W boson production in heavy ion collisions might lead to an improved determination of 
the nPDEs. Moreover, asymmetries in the individual yields of and should permit the 
flavour decomposition of u and d quark distributions in nuclei. 

The ATLAS ||8l 0 and CMS HOl |TT1 collaborations reported the observations of Z bosons in 
heavy ion interactions, at a centre-of-mass energy of 2.76 TeV per nucleon pair. These data 
showed that the Z boson yields per nucleon-nucleon (NN) collision are essentially unmodified 
by the medium produced in the collisions. Although W bosons decaying to a lepton and a 
neutrino are more difficult to detect, their rate is about ten times larger than that of Z bosons 
decaying to leptonic final states. The production of W bosons in PbPb collisions was reported 
by CMS HT^ and ATLAS Ifl^ . using data corresponding to an integrated luminosity of 7.3 
and 150 collected in 2010 and 2011, respectively. The W boson yield per NN collision 
was shown to be compatible with the one measured in pp collisions, when taking into account 
isospin effects arising from the mixture of protons and neutrons in the colliding nuclei. How¬ 
ever, the presence of 10-20% nPDL effects on Z and W boson production could not be excluded 
due to the relatively large experimental and theoretical uncertainties of these results. 

The 2013 pPb LHC run provides the best currently available data sample to look for initial-state 
effects (such as PDE modifications) using electroweak bosons. The NN-equivalent luminosity 
is of the same order of magnitude as for the 2011 PbPb run, and the production cross sections 
are approximately a factor of two greater owing to the increased energy, 5.02 TeV per nucleon 
pair. Lurthermore, the asymmetry of the pPb collision system allows for the measurement of 
other observables such as forward-backward pseudorapidity asymmetries. This Letter reports 
a study of W boson production in a sample of pPb collisions corresponding to an integrated 
luminosity of (34.6 ± 1.2) nb^^ HH, collected by the CMS experiment. 


2 Experimental methods 

The direction of the proton beam was initially opposite to the positive direction of the CMS 
longitudinal axis |[l5l , and was reversed after 60% of the data were taken. The beam energies 
were 4 TeV for protons and 1.58 TeV per nucleon for lead nuclei, resulting in a centre-of-mass 
energy per nucleon pair of = 5.02 TeV. As a result of the energy difference of the colliding 
beams, the NN centre-of-mass frame in pPb collisions was not at rest with respect to the labora¬ 
tory frame. Massless particles emitted at pseudorapidity rj in the NN centre-of-mass frame are 
detected at ?/iab = V ~ 0-465 (first proton beam orientation) and J/iab = ?/ + 0.465 (second proton 
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2 Experimental methods 


beam orientation) in the CMS coordinate system, as defined in Ref. llTSl . The results presented 
hereafter are expressed in the usual convention of the proton-going side defining the positive 
pseudorapidity It coincides with the CMS convention in the second period of data taking, the 
first one being reversed before summing yields from the two beam configurations. 

A detailed description of the CMS detector can be found elsewhere (TSl. Its central feature 
is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. 
Within the field volume are the silicon pixel-and-strip tracker, a lead tungstate crystal electro¬ 
magnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL), each 
composed of a barrel and two endcap sections. The silicon tracker consists of 66 M pixel and 
10 M strip sensor elements, and measures charged-particle trajectories in the pseudorapidity 
range |f/iab| < 2.5. Outside of the solenoid, muons are detected in the |?/iab| < 2.4 range, with 
gas-ionization detector planes based on three technologies: drift tubes, cathode strip cham¬ 
bers, and resistive-plate chambers. Electrons are identified in the ECAL, which is made of 
75 848 lead tungstate crystals and covers \tj\ah\ < 1-48 in the barrel and 1.48 < \tj\ah\ < 3.00 
in the two endcap regions. The CMS apparatus also has extensive forward calorimetry, in¬ 
cluding two steel/quartz-fiber Cherenkov hadron forward (HE) calorimeters, which cover the 
2-9 < |?/iab| < 5.2 range. Eor online event selection, CMS uses a two-level trigger system. 

Selection criteria similar to the ones developed in Ref. IT^ are applied to the pPb sample to 
remove events with electromagnetic, beam-gas, or multiple collisions (pileup). The W boson 
yields are corrected for the induced (4.0 ± 0.5)% signal loss. 

The primary signature of a W boson is a high transverse momentum (px) lepton. The current 
analysis is restricted to leptons of px greater than 25GeV/c. The muon analysis is based on a 
sample triggered by requiring a single muon with px above 12 GeV/c, while the electron ana¬ 
lysis uses an ECAL-triggered sample with a transverse energy threshold of 15 GeV. Leptons 
are reconstructed with the same algorithms as in proton-proton collisions lllTl[l8ll . and stan¬ 
dard selection criteria are applied, as in Refs. |[l2l[T9l. A special electron charge determination, 
as described in Ref. Il20l , is used in order to reduce the electron charge misidentification to a 
sub-percent level. Events are reconstructed using particle-flow (PE) techniques II2T11221 . which 
reconstruct and classify individual particles with an optimised combination of all subdetector 
information. 

Two criteria are used to remove specific background sources. Eirst, events with two oppositely 
charged leptons, with the second lepton px greater than 15 (10) GeV/c for muons (electrons) are 
removed, since they correspond to well-identified processes like Drell-Yan, Z boson or high-px 
quarkonium production. Second, the leptons are required to be isolated, in order to reduce the 
contamination coming from jet fragmentation. The energies of all PE candidates are summed 
within a cone centred around the lepton, with the exception of the lepton itself. The lepton is 
considered isolated if the total transverse energy in the cone is small compared to its transverse 
momentum. Eor muons, a cone of radius AR = \/ (Ap)^ -|- {A(p)^ = 0.3 is used, where At] and 
A(p are the pseudorapidity and azimuthal distances to the lepton. The candidate is rejected 
if the in-cone transverse energy is greater than 10% of the muon px. Eor electrons, a cone 
of AR = 0.4 is used, and only particles with px greater than 1 GeV/c are summed, to reduce 
the underlying-event enhanced contribution. The electron candidate is rejected if the resulting 
transverse energy is greater than 11.5% (9.5%) of the electron px, for the EGAL barrel (endcaps). 

An important characteristic of events containing a W —)> £v decay is the missing transverse 
energy (Ex) associated with the undetected neutrino. It is computed as the magnitude of the 
vectorial sum of transverse momenta of all the PE candidates in the event. The analysis is 
performed using ten lepton pseudorapidity bins, each 0.5 wide except for the most forward 
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and backward regions (2 < \rj\a\)\ < 2.4). After having allied the lepton selection criteria, 
examples of the resulting distributions are shown in Fig.^for and e+, in the most central 
(—0.5 < J/iab < 0.0) and furthest forward (2.0 < ;^iab < 2.4) ranges. The distributions for other 
bins and for the negative leptons are similar. 
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Figure 1: Missing transverse energy distribution for W+ —> }i^v (left) and W+ —?■ e+r (right) 
events within the —0.5 < J/iab < 0.0 (top) and 2.0 < //lab < 2.4 (bottom) ranges. Binned fits 
to the data (red points) are performed with four contributions, stacked from bottom to top: 
multijet (QCD, blue), W+ —> r+r (brown), Z —> (white) and W+ —)■ t^v (yellow). The 

regions are defined such that the proton is moving towards positive r]\ah values. Error bars 
represent statistical uncertainties. The lower panels display the data divided by the result of 
the fit, with the band representing the statishcal uncertainties on the sum of the fit components, 
for each bin. 


To extract the number of events with a lepton coming from a W boson, binned fits of these 
distributions are performed, including the signal and main background contributions, in each 
riiab bin. The shapes assumed for the electroweak processes, namely the —)• £^v signal 

as well as background from —> r^r and Z —)■ are determined by the simulations 

described hereafter, taking into account the acceptance and efficiency. Their relative normal¬ 
ization is given by the unmodified theoretical cross sections (as computed in Ref. ||23j). A max¬ 
imal 20% variation of the W/Z normalization ratio is taken into account, due to potentially 
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2 Experimental methods 


different nuclear modifications of the Z and W bosons, and resulting in a 1-3% systematic un¬ 
certainty in the extracted W yields. The noticeable difference between the distributions for 
the Z —> e+e^ and Z —> processes in the forward region (bottom plots of Fig.[^ results 

from the greater ECAL coverage allowing missed electrons with 2.4 < \rj\ah\ < 3.0 to be ac¬ 
counted for in the Ex calculation. The shape of the QCD multijet background is modelled by 
the functional form/(Ex) = (Ex + Ex'^)“ + Et°)- It is shown to reproduce the Ex 

shape of data events containing non-isolated leptons, with the Ex°/ and /3 parameters, which 
are observed to depend mildly and linearly on the cone/lepton transverse energy ratio. These 
fitted parameters are then extrapolated to the isolated lepton signal regime and the resulting 
function is used as the QCD background shape. The multijet background contribution is larger 
in the electron channel because the misidentified lepton rate is higher, particularly due to a 
contribution from photon-jet events. Contributions from other sources, such as tt production 
and high-px quarkonia, were found to be negligible. 

A small charge misidentification correction (less than 0.2%) is applied to the electron yields; 
this correction is negligible for muons. All fits are of good quality, as illustrated by the bottom 
panels of Fig. that show the ratio of the data to the fit outcome. The observed numbers of 
leptons coming from W boson decays over the entire pseudorapidity range are: 11660 ±111 
9459 ± 99 ]i^, 9892 ± 116 e+, and 7872 ± 101 or, where the uncertainty is statistical, determined 
by the fit procedure. 

In order to correct for inefficiencies in the lepton trigger, reconstruction, and selection, the elec- 
troweak processes W —> £v have been simulated using the pythia 6.424 generator |24^1 with 
a mixture of pp and pn interactions corresponding to pPb collisions. The detector response to 
each PYXHIA signal event is simulated with Geant4 ESI and then embedded in a minimum 
bias pPb background event. These background events are produced with the HIJING event 
generator I26l and passed through Geant4 as well. Each simulation is done twice, once for 
each proton beam direction, and includes a boost to reproduce the 0.465 rapidity shift. The 
embedding is done at the level of detector hits, and the signal and background events share 
the same generated vertex location. The embedded event is then processed through the trigger 
emulation and the full event reconstruction chain. The resulting reconstructed events are then 
reweighted to match the distributions observed in data of the event vertex and activity (as mea¬ 
sured in the HE calorimeters). The obtained efficiencies vary with (with higher efficiencies 
at mid-rapidity), from 59% to 89% for muons, and from 51% to 84% for electrons. 

The various components of the single-lepton efficiency are also directly computed from pPb 
data, using 7. ^ It samples, and techniques described in Ref. II23H . These efficiencies are then 
compared to the corresponding efficiencies computed from simulations. In the case of trigger 
and reconstruction efficiencies, they are found to be consistent. The isolation criterion rejects 
more leptons in data, because the local activity of the underlying event is greater than in the 
simulation. To account for such discrepancies, the efficiency from W —> simulation is mul¬ 
tiplied by correction factors, which are determined as the ratio of the single-lepton efficiencies 
measured in Z —> data to those estimated in simulations. The so-called "tag-and-probe" 
method used for this estimation is described in Ref. |[27]. These correction factors are computed 
in bins of //lab and for positively and negatively charged muons separately. In the electron case, 
the low statistical precision motivates a correction factor estimated for electrons and positrons 
combined. 

The total systematic uncertainty in the lepton yields is estimated by adding the different con¬ 
tributions in quadrature. The /^lab-dependent sources of systematic uncertainty arise from the 
method used for the estimation of multijet background (0.1-2.0% for muons, 0.5-3.8% for elec- 




5 


trons), the normalization of the electroweak background (1-3% for muons and electrons), the 
efficiency correction factors (2.2-7.5% for muons, 2.6-7.4% for electrons), and the energy scale 
of electrons (0.1-2.0%). The uncertainty in the momentum scale of muons is found to be negli¬ 
gible. The integrated luminosity measurement uncertainty (3.5% 1141 ) affects only the W boson 
production cross sections and cancels in the asymmetry measurements, as does the additional 
global uncertainty arising from the efficiency of the filter rejecting pileup events (0.5% for both 
channels). Though the common electron/positron correction factors cancel, a residual system¬ 
atic uncertainty of 3% is assigned to the charge asymmetry, based on simulation studies and 
j/iab-integrated efficiencies determined from Z —e+e^ data. No other systematic uncertainty 
cancellations are assumed for the asymmetry results. 


3 Results 

Figure l^shows the production cross sections for pPb —> -|- X —)■ -|- X as a function of the 

charged lepton pseudorapidity in the laboratory frame, with the lepton having pj > 25 GeV/c. 
The cross sections are determined by dividing the efficiency-corrected lepton yields by the 
integrated luminosity. 



^lab ^lab 


Figure 2: Production cross sections for W+ —)■ (left) and —>• i~v (right), as a function 

of the lepton pseudorapidity. Error bars represent the statistical uncertainties, while brackets 
show statistical and systematic uncertainties summed in quadrature. The global luminosity 
uncertainty of ±3.5% is not included. To improve visibility, the muon (electron) measurements, 
in red circles (blue squares), have been shifted by —0.05 (±0.05) in pseudorapidity. 

Since the cross sections measured in the electron and muon channels are found to be in good 
agreement with each other, they are combined using the BLUE method [ 281 . Figure|^compares 
the combined cross sections with next-to-leading-order (NLO) perturbative QCD predictions 
provided by the authors of Ref. Hi using CTIO Il29l proton parton distribution functions (PDF) 
without or with EPS09 Il30ll nPDF corrections, termed CTIO and CT10+EPS09, respectively. 
Their uncertainties are estimated as prescribed in Refs. l30l . Table gives the measured 
cross sections for each channel separately and combined, as a function of the lepton pseudo¬ 
rapidity, for positive and negative leptons. The theoretical predictions and their uncertainties 
(coming from the PDF set and from the renormalisation and factorisation scales) are also given. 
The agreement between the data and both theoretical predictions is within the uncertainties, 
although a small excess of candidates appears at negative J/iab, i-e. in the Pb ion beam 
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3 Results 


direction. 

The comparison between the CTIO and CT10+EPS09 calculations shows that the predicted 
modifications of the PDFs are of the same order as the theoretical uncertainties. This indi¬ 
cates that cross sections alone lack discriminating power, and motivates the study of various 
asymmetries of the and cross sections. The interest in such asymmetries is twofold. First, 
some of the experimental (e.g. integrated luminosity) and theoretical (e.g. scale dependence) 
uncertainties cancel in such asymmetries. Second, the various asymmetries exhibit different 
sensitivities to the nuclear modifications of the PDFs, as discussed below. 
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Figure 3: Production cross sections for (left) and W~ —>• v (right), as a function 

of the lepton pseudorapidity. Error bars represent the statistical uncertainties, while brackets 
show statistical and systematic uncertainties summed in quadrature. The global luminosity un¬ 
certainty of ±3.5% is not displayed. Theoretical predictions with (CT10+EPS09, dashed green 
line) and without (CTIO, solid red line) PDF nuclear modifications are also shown, with the un¬ 
certainty bands. The bottom panels show the ratio of the data (black points) and CT10+EPS09 
(dashed green line) to the CTIO baseline. All theory uncertainty bands include scale and PDF 
uncertainties, except the EPS09 of the bottom panels which only includes the EPS09 PDF un¬ 
certainties. 


The lepton charge asymmetry, defined as (N^ — ) / (N^ ± ) with NJ*" being the efficiency- 

corrected lepton yields, is shown in Fig. as a function of J/iab, and compared to the theoret¬ 
ical predictions. For > ~1/ both calculations reproduce the present measurements. For 
'/lab < —1/ however, the two calculations overpredict the asymmetry values. A possible physi¬ 
cal origin of this disagreement could be a different modification of u and d quark distributions 
in nuclei. In proton-(anti)proton collisions, the W-boson charge asymmetry is known to be a 
sensitive probe of the down-to-up quark PDF ratio in a proton, d^/uP 120113111321. Similarly, 
this asymmetry in pPb collisions measured in the lead fragmentation region (i.e. J/iab < 0.465) 
probes these quark densities in a nucleon inside the lead nucleus. Assuming the standard 
isospin symmetry (uP = d", u” = dP), one can define a similar ratio, = dP/uP x 

Rd/^u/ where R; are the nPDF ratios, Ru = uP^^/uP and R^ = dP^^/d^. The typical quark 
momentum fraction probed in the Pb nucleus is given by x ~ Mw/x exp(—j/jab + 0.465) 
(assuming that the W boson rapidity is similar to that of the lepton), therefore x ~ 0.02 — 0.20 
in the range —2 < rji^h < 0- In most global fit analyses of the nPDFs (as in the case of EPS09), 
it is assumed that the nuclear ratios respect the isospin symmetry, namely Ru = Rd, essentially 
to minimise the number of free parameters in the fits. However, no physical reason prevents 
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Table 1: Production cross section for pPb —>W + X—)-£i/ + X for positively (top) and nega¬ 
tively (bottom) charged leptons of pj larger than 25 GeV/c, in nanobarns, as a function of the 
lepton pseudorapidity Values are given first for muons and electrons separately, then com¬ 
bined. Quoted uncertainties are first statistical, then systematic. Theoretical predictions with 
(CT10+EPS09) and without (CTIO) PDF nuclear modifications are also given, with their uncer¬ 
tainties. The global normalization uncertainty of 3.5% is not included in the listed uncertainties. 


If (nb) [^bin] 

[-2.4,-2.0] 

[-2.0,-1.5] 

[-1.5,-1.0] 

[-1.0,-0.5] 

[-0.5,0] 


43.0 ±2.2 ±3.1 

62.5 ± 2.1 ± 2.6 

86.9 ±2.6 ±3.4 

98.1 ±2.7 ±2.6 

98.3 ±2.8 ±3.3 

e+ 

46.5 ± 2.6 ± 3.6 

64.0 ±3.1 ±4.2 

84.2 ±3.1 ±4.8 

99.8 ±3.0 ±4.6 

102.0 ±2.9 ±4.6 

e+ 

44.5 ±1.7 ±2.3 

62.9 ±1.8 ±2.2 

85.9 ±2.0 ±2.7 

98.6 ±2.1 ±2.3 

99.7 ±2.1 ±2.7 

CT10+EPS09 

42 1+2-6 

66 . 0 ++® 

84.6-;® 

9 X 4+53 

96.01®;® 

CTIO 

43.4+1® 

65.8+®;^ 

82.4—;® 

90.511® 

94.4l®;J 

If (nb) Dbin] 

[0,0.5] 

[0.5,1.0] 

[1.0,1.5] 

[1.5,2.0] 

[2.0,2.4] 


113.9 ±3.1 ±4.5 

101.3 ±2.8 ±2.9 

102.3 ± 2.8 ± 3.6 

107.9 ±3.1 ±5.7 

107.8 ±3.7 ±8.4 

e+ 

99.6 ±2.7 ±3.6 

102.8 ±2.9 ±4.6 

95.6 ±3.4 ±5.8 

95.4 ± 3.5 ± 6.2 

108.3 ± 4.3 ± 8.7 

i+ 

105.3 ±2.1 ±2.8 

101.8 ±2.1 ±2.5 

100.2 ±2.2 ±3.1 

102.3 ±2.3 ±4.2 

108.1 ± 2.8 ± 6.0 

CT10+EPS09 

QC 0 + 6.2 

^^•^-6.4 

95(5+11 

95(7™ 

95(3-® 

91(61- 

CTIO 

97.0+|8 

100.0l|4 

103.411® 

105.7™ 

103.61®;° 

% (nb) [rj bin] 

[-2.4,-2.0] 

[-2.0,-1.5] 

[-1.5,-1.0] 

[-1.0,-0.5] 

[-0.5,0] 


74.5 ± 3.0 ± 5.6 

84.5 ± 2.8 ± 4.4 

89.4 ±2.6 ±3.5 

81.4 ±2.5 ±2.6 

80.6 ± 2.6 ± 2.6 

e“ 

70.2 ± 3.2 ± 4.8 

74.3 ± 3.3 ± 4.8 

79.6 ± 3.1 ± 4.3 

80.7 ±2.7 ±3.7 

81.3 ±2.6 ±4.0 

t- 

72.1 ± 2.2 ± 3.7 

79.9 ± 2.1 ± 3.3 

85.4 ±2.0 ±2.7 

81.1 ± 1.8 ± 2.1 

80.8 ±1.9 ±2.2 

CT10+EPS09 

65.2-- 

72.4^- 

75.9—® 

76(9-® 

76.11- 

CTIO 

64.2+3;9 

70.1-:7 

73 3+^-^ 

/d.d_4 8 

74.81- 

75.11^;^ 

If (nb) [rj bin] 

[0,0.5] 

[0.5,1.0] 

[1.0,1.5] 

[1.5,2.0] 

[2.0,2.4] 


81.7 ±2.5 ±3.0 

78.8 ± 2.5 ± 3.3 

69.8 ±2.3 ±3.0 

62.9 ±2.1 ±3.3 

63.1 ±2.8 ±5.1 

e“ 

73.5 ± 2.5 ± 3.5 

74.0 ± 2.5 ± 3.5 

70.6 ±2.8 ±4.6 

55.0 ±2.7 ±4.1 

64.6 ± 3.3 ± 6.0 

t- 

78.0 ±1.8 ±2.3 

76.5 ± 1.8 ± 2.4 

70.1 ±1.8 ±2.5 

59.8 ±1.7 ±2.6 

63.7 ±2.1 ±3.9 

CT10+EPS09 

73.6^1+ 

69.7-- 

64.8-® 

5 V 1 -® 

53 4+^-^ 

CTIO 

74.3^« 

72.4—;® 

69.1-® 

64.51®;® 

59 3+^-^ 


nuclear modifications to be different for up and down quark PDFs. For example, it is known 
that the shapes of the up and down quark distributions in protons are different (331. Further¬ 
more, the present disparity between data and theory is unlikely to come from the proton PDF 
assumption, given the excellent agreement of lepton charge asymmetry measured in pp colli¬ 
sions by CMS (321 and ATFAS (Ml with NFO calculations using CTIO parton densities. 

A traditional way to probe nuclear parton densities is to compare the forward and backward W 
yields, that are respectively sensitive to the nPDFs at small and large x. The forward-backward 
asymmetries {+p\ah)/ (—?/iab) are shown in Fig.separately for the positively and neg¬ 
atively charged leptons, and compared to the same predictions as mentioned above. Given the 
experimental accuracy and the magnitude of the differences between the two sets of predic¬ 
tions, the measurements have a potential to discriminate between them. However, although 
the negative lepton decay channel appears to slightly favour the CT10+FPS09 prediction over 
the CTIO calculation, the positive lepton channel does not, thus no firm conclusion can be 
drawn. 

Another asymmetry variable, {N^{+r]iab) - N/(-?7iab))/(N‘7(-F//iab) - ^^(-T/iab)), was pro¬ 
posed in Ref. (H to reach maximum sensitivity to nuclear modifications of PDFs. However, 
this asymmetry probability distribution shows a very non-Gaussian behaviour, when its de¬ 
nominator approaches zero, and its sign can be flipped within the uncertainty. A different 
asymmetry is proposed here, Ng{+ri\ah)/N£{—ri\ah), a forward-backward asymmetry of the 
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3 Results 
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Figure 4: Lepton charge asymmetry, (N^ — ) / (N^ + ), as a function of the lepton pseu¬ 

dorapidity Error bars represent the statistical uncertainties, while brackets show statistical and 
systematic uncertainties summed in quadrature. Theoretical predictions with (CT10+EPS09, 
dashed green line) and without (CTIO, solid red line) PDF nuclear modifications are also 
shown, with their uncertainty bands. 
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Figure 5: Forward-backward asymmetries, N£(-|-)^iab)/N^(—?/iab)/ for the positive (left) and 
negative (right) leptons. Error bars represent the statistical uncertainties, while brackets show 
statistical and systematic uncertainties summed in quadrature. Theoretical predictions with 
(CT10+EPS09, dashed green line) and without (CTIO, solid red line) PDF nuclear modifications 
are also shown, with their uncertainty bands. 
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charge-summed W bosons, which achieves a similar sensitivity. As in the case of the charge 
asymmetry, this asymmetry can be related to the nuclear modifications of the PDFs within the 
lead nucleus. Here, forward (backward) W boson production is sensitive to the PDFs of the sea 
quark at x ~ 10^^ (valence quark at x ~ 10^^) in the lead nucleus. Therefore, the forward- 
backward ratio probes the small-x modification of the lead nucleus PDF (shadowing) over the 
large-x modifications (anti-shadowing). This asymmetry is shown in Fig. and deviates from 
unmodified PDFs, more clearly favouring CT10+EPS09 over CTIO. 



Figure 6: The forward-backward asymmetry of charge-summed W bosons, as a function of the 
lepton pseudorapidity. Error bars represent the statistical uncertainties, while brackets show 
statistical and systematic uncertainties summed in quadrature. Theoretical predictions with 
(CT10+EPS09, dashed green line) and without (CTIO, solid red line) PDF nuclear modifications 
are also shown, with their uncertainty bands. 

In order to quantify the agreement between the data and the expectation from the CTIO and 
CT10+EPS09 calculations, a ^ test is performed for each of the above (correlated) variables. 
The few correlations in experimental uncertainties described above, only relevant for boson 
cross sections but not for asymmetries, are taken into account, as well as the correlations in 
theoretical uncertainties. The resulting values and probabilities are given in Table The 
CT10+EPS09 calculations provide a better description of the data, with still a relatively low 
probability for the lepton charge asymmetry, because of the backward region. 


4 Summary 

The first measurement of W boson production in pPb collisions has been reported, using the 
electron and muon decay modes for leptons of px above 25GeV/c and |piab| < 2.4. The dif¬ 
ferential cross sections as a function of the lepton pseudorapidity agree with theoretical pre¬ 
dictions assuming both unmodified (CTIO) and modified (CT10+EPS09) nPDFs, except in the 
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4 Summary 


Table 2: Values of the test between the measurements and the theoretical predictions, with 
(CT10+EPS09) or without (CTIO) nuclear modifications of the PDFs. The probability (Prob.) to 
measure a value greater to that measured in data is also given for ten degrees of freedom in the 
case of the first three observables and five degrees of freedom for the three others observables. 


Observable 


GTIO 

GT10+EPS09 


Prob.(%) 


Prob.(%) 

d(7/d//(f'+) 

13 

25 

8.6 

57 

da/dr]{£^) 

15 

14 

8.2 

60 

(N+-N,-)/(N++N,-) 

15 

12 

11 

35 

(“l“?/lab) / { ^lab) 

3.1 

68 

3.2 

68 

Nf (+?/lab)/Nf (-?/lab) 

9.7 

8.4 

3.5 

63 

^^ei+Vlab) /Ne{ — Vlab) 

6.2 

29 

2.1 

83 


most backward region (Pb ion beam direction), where a hint of an enhancement is seen for the 
bosons. In the same region, the related lepton charge asymmetry deviates slightly from 
the predictions, something that could potentially arise from different nuclear modifications of 
the up and down quark PDFs. In a related observation, forward-backward asymmetries show 
a deviation from unmodified PDFs. Taken together, these measurements show the need for 
including W boson data in nuclear parton distribution global fits. 
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